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D 

pH 1 We construct commuting transfer matrices for models describing the interaction between a single 

quantum spin and a single bosonic mode using the quantum inverse scattering framework. The trans- 
fer matrices are obtained from certain inhomogeneous rational vertex models combining bosonic and 
spin representations of SU(2), subject to non-diagonal toroidal and open boundary conditions. On- 
ly open boundary conditions are found to lead to integrable Hamiltonians combining both rotating 
and counter-rotating terms in the interaction. If the boundary matrices can be brought to trian- 



O 
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gular form simultaneously, the spectrum of the model can be obtained by means of the algebraic 
Bethe ansatz after a suitable gauge transformation; the corresponding Hamiltonians are found to 

rSJ , be non-hermitian. Alternatively, a certain quasi-classical limit of the transfer matrix is considered 

> ■ 

where hermitian Hamiltonians are obtained as members of a family of commuting operators; their 
CO , diagonalization, however, remains an unsolved problem. 

CO 

o 

|^ ' I. INTRODUCTION 

o ' 

-I— > 

Models describing the interaction between bosonic modes and spin degrees of freedom play an important role 
in physics. The prototype application is provided by atom-radiation interactions 1], where the atomic dipole is 
described by an effective spin interacting with a single mode electric field . More recently, this class of models 

& ■ 

has been considered to investigate various aspects in mesoscopic physics where a given two-level system interacts 

o ■ n 

• • . with a bosonic environment 41. Furthermore, we mention applications to analyze the decoherence in superconducting 

... . Q 

circuits [5|, |6J, and to describe ions in harmonic traps [7]. Both these quantum devices are of potential importance for 
quantum computation Finally, applications have emerged in quasi-2D semiconductors in transverse magnetic 

— ' field HQQ. 

In this paper we focus on models describing a single spin S interacting with a single mode boson field which were 
introduced originally by Jaynes, Tavis and Cummings P, In two cases the corresponding model Hamiltonian is 



known to be exactly diagonalizable with elementary means: for S — *■ oo where the model is quadratic 13H and for weak 

nnn 

or resonant interactions, where the Rotating Wave Approximation (RWA) can be applied [2, 3, 14]. In the traditional 
framework of laser physics, the RWA assumes that the relevant dynamics is given by coherent oscillations of the 
population of the atomic energy levels. Therefore only photon emissions accompanied by an excitation of the atom 
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(and vice versa), described by operators of the type S~ , are taken into account in the atomic dipole coupling with 
electric field E ■ d ~ (a + a))S x . The so-called counter- rotating terms a< S + + aS~ , although contained in the coupling, 
can be neglected at resonant cavity frequency Qo since they induce a short time dynamics (compared to the rotating 
terms aS + + a^S^); as a result of the approximation the operator S z + a^a is conserved. Such a picture is robust 
for weak fields and as long as the cavity mode can be tuned conveniently. In the most recent applications, however, 
one or even both these conditions are not met. In superconducting circuits, for example, the effective spin-boson 
coupling to the electromagnetic field of the device can be naturally large, thereby involving multi-boson processes, 
which are characteristic for the presence of counter rotating terms 15[ ; furthermore numerical inspections indicate 
that such terms can bend the energy levels whose slope should be controlled reliably to study how the noise affects 
the coherence in the system Finally we mention that in semiconductors, rotating and counter rotating terms 
result from Rashba and Dresselhaus spin-orbit interactions; the physical origin of the spin transport in such systems 
is an open problem that is intensively studied 12 1 . 



For finite spin representations and generic spin-boson interaction the Hamiltonian is a non trivial matrix of infinite 
rank. The main complication comes from the absence of a conserved "number operator" of the type S z + a^a as in the 
RWA that would have made the problem finite dimensional - a^a is the bosonic number operator. In this sense, using 
the language of spin chains, the problem will have the complexity of an XYZ chain. We will speak of a spin-boson 
model to contain rotating terms only, if the Hamiltonian action only connects a subspace of finite dimension. It is 
worth noting that it is the simultaneous presence of so called rotating, S^a, and counter-rotating terms S'a' that 
makes the Hamiltonian non-trivial in the sense that a reduction to a finite dimensional subspace is impossible. 

In this work we study integrable spin-boson models that can be constructed by means of the quantum inverse scat- 
tering method (QISM) [17| . The central object is the transfer matrix, that it is a generating function of integrals of the 
motion, including the Hamiltonian by definition. The formalism includes a certain freedom to choose (twisted) peri- 
odic boundary conditions and open boundary conditions, where two boundary matrices can be included, representing 
boundary magnetic fields in the standard formalism for spin chains [Lsl [l^ | . Twisted boundary conditions have been 



investigated for XXX-symmetry, i.e. for models emerging from the XXX i?-matrix in [221 ] and hermitian spin boson 
models have been created that still have a conserved number operator, indicating that the apparent counter-rotating 
terms are spurious. In the present work, open boundary conditions for XXX-type symmetry are analyzed, where 
two non-parallel boundary matrices are meant to break the XXX symmetry completely. Integrable open boundaries, 
corresponding t o dy namic boundary matrices (i.e. depending on the spectral parameter) obey the reflection equation 
due to Sklyanin 19(. Two approaches are used. First, a bosonic Lax matrix, relative to the XXX i?-matrix is employed 
and the QISM and algebraic Bethe ansatz is elaborated for this case. Second, the spin-boson Hamiltonian is obtained 
from a two-site spin chain with mixed representations by virtue of a 'large S limit' (algebraic contraction) of one of 
the spins to a bosonic degree of freedom. 

The paper is laid out as follows. In the next section we summarize the ingredients of the reflection equations to 
deal with integrable models with generic boundaries. In the Section IIIII we present the spin boson transfer matrix. 
The diagonalization will be pursued by a 'direct' algebraic Bethe ansatz (Section llVAp and by means of the algebraic 
contraction of a suitable auxiliary spin problem (Section HVB[) . The Hamiltonian models extracted from the transfer 
matrix will be presented in the Section [Vj The Section PVTl is devoted to the conclusions. 
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II. INTEGRABLE BOUNDARY CONDITIONS 

In this section, we review the basic aspects of the quantum inverse scattering method providing integrable theories 
with non-trivial boundary conditions. 

The central object of the method is the transfer matrix i(A), which can be conveniently written as the trace over 
an auxiliary space A of an ordered product of L distinct Lax matrices Cj(X) acting in the site j. It can be explicitly 
written as 

t(X) = tr A [T(\)}, T(X)=jC l (X)jC l . 1 (X)...jC 1 (X), (1) 

where T(A) is called monodromy matrix. 

The i?-matrix and Lax matrices Cj(X) are solutions of the Yang-Baxter equations 

R 12 (X-X')C ( ^(X)Cf\x') = cf\x')Cf\x)R 12 {X-X'), (2) 
Ri2(X)Ri 3 (X + X')R 23 (X') = R 23 (X')R 13 {X + X')R 12 (X) , (3) 

where the suffix m in indicates that in a multiple tensor product of auxiliary spaces, the Lax matrix acts 

non-trivially only on the m-th copy. A rational solution of ([3]) is the following R matrix 

/ X + r/ \ 

A 77 

R(X)= , (4) 

77 A 

\ A + 77 / 

where 77 is a C — number named quasi-classical parameter. The Yang-Baxter equation for the Lax matrices © provides 
the condition to the commutation property of the transfer matrix (fT]) and give rise to the fundamental commutation 
relations (FCR) for their operator entries. 

In the realm of the QISM boundary terms can be considered. The simplest way to do it is to construct the 
monodromy matrix in terms of modified Lax matrices C(X) — > QC{X), Q being a matrix of constants (with respect to 
the spectral parameter) [18j|. The corresponding theory is still integrable so far the quadratic relation §Z§ is satisfied. 
For that it suffices that 

[#i 2 (A),£ (1) £ (2) ]=0 (5) 

For the case of i?-matrix (T4|) this relation is satisfied by any 2x2 matrix with entries g\\, g± 2 , g 2 \, g 22 independent of 
the spectral parameter. 

It turns out that such modified monodromy matrices correspond to theories with toroidal boundary conditions. 
In the case where the boundary matrix Q is diagonal the corresponding transfer matrix (fT]) can be diagonalized 
with very little difference from the periodic case because it does not change in a drastic way the properties of the 
monodromy matrix elements. In general this does not occur when Q is non-diagonal. Nevertheless, for the rational 
solutions of the Yang-Baxter equation (J3j) , it is possible to apply a combination of suitable auxiliary and quantum space 
transformations which relate the non-diagonal boundary problem with a diagonal one [27] . These transformations are 
the irreducible representations of the symmetry of the auxiliary and quantum space of the ^-operator Cj (A) . Such 



invariance of the Lax matrix can be expressed by [2 



4° (A) 



4 13 (A) 



gwu< 



(6) 



where Uj stand for the irreducible representation of the symmetry of the quantum space 27j. For instance, if the 
quantum space is invariant by spin-S" representation of SU(2), the relation © is satisfied by 



us,s us,s-i ■ ■ ■ us-s 



V U-S,S U-S,S-1 



(7) 



-s,-s J 



where some of the matrix elements 



2Sj-k k 
US,S-k = 9l l { 9l,2\ 



(25)! 



k\{2S-k)V 



(8) 



while the remaining satisfy the following recurrence relation 

(51,152,2 - 31,232,1) 



US-k-l,S-r 



71(25-71+1) 




fc(25-*+l) 



MS-fc+l.S- 



y (fe + i)(2S-fe) 

32,1 2(S-k) 
3i,i v/(fc + l)(25-fe) 



"S-fc,S-n+l 



US-k,S- 



(9) 



(k + l)(2£-fc) 

and /c, n are integers satisfying fc,n = 0,...,25. 

As the transfer matrix ([TJ) gives origin to a commutative family of operators, it can be regarded as generating 
functional of the integrals of the motion of the quantum theory, among which one of them can be taken as the 
Hamiltonian. For example the logarithmic derivative of i(A) provided a protocol to demonstrate the integrability of 
a variety of local models (that is with nearest neighbor interactions) 17] . For later use we remark that an arbitrary 
function of t(A) or e.g. its coefficients in a Taylor series in A can be taken to be the Hamiltonian. 

A different class of lattice models can be obtained focusing on the quasi-classical parameter 77. Specifically, integrable 
long-range Hamiltonians, as the Gaudin and the BCS model, emerge from the so called quasi classical limit 77 — * of 
vertex models 123. 



24 



25 
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More complicated boundary terms can be considered in the integrable theory via a procedure demonstrated by 
SklyaninQQ. To begin with, let's impose the following constrains on the i?-matrix 



Unitarity: fli 2 (A)i?i 2 (-A) = C(A)1; 
Parity invariance: Pi2-Ri2(A)Pi2 = R12W; 
Temporal invariance: i?i2(A)* lt2 = i?i2(A); 
Crossing symmetry: i?* 1 2 (A) J R* 1 1 2 (-A - 2p) = C(A + p)l; 



(10) 

(11) 
(12) 

(13) 



where C(A) = t? 2 — A 2 and the crossing parameter p = r/. Here 11 is the identity matrix, P12 is the permutation operator, 
t a denotes transposition on the a-th space. 

According to the Sklyanin procedure, open boundary conditions preserve the integrability of the model need to 
provide 2-particle scattering in the bulk 'compatible' with the scattering off the boundaries. Such boundary conditions 
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can be parametrized by boundary matrices K(X) which satisfy the so called reflection equations 

Ri 2 (X - A') K- (A)i? 21 (A + A') K- (A') =K- (X')R 12 (X + A') K- (X)R 21 (X - A'), 

( 14 ) 

R 21 (-X + A') K\ (X)R 12 (-X - A' - 2rj) K\ (A') =K\ (X')R 21 (-X - A' - 2 V ) K\ (X)R 12 (-X + A'), 

1 2 

where K (A) = K(X) ® 1, if (A) = 1 ® ^(A), if is the transpose of if, and i?i2 is a solution of the Yang-Baxter 
equation ([3]). 



General non-diagonal C-number solutions in case of the rational i?-matrix ^ are the if-matrices 31 1 




K±(X) = • (15) 



Within this framework the analogue of the row-to-row transfer matrix ([T]) as generator of commuting integrals of 
motion in the case of open boundaries is the following operator 

f(A)=tM[iC+(A + ry)[T(-A)]- 1 K_(A)T(A)] . (16) 

often called 'double row transfer matrix' in the literature. For diagonal AT-matrices the diagonalization of t(X) can 
be done by an extension of the algebraic Bethe Ansatz [ijj]. For models based on the rational i?,-matrices (@| con- 
sidered in this paper this approach can be extended to certain non-diagonal if-matrices by combination of suitable 



transformations in both quantum and auxiliary space [3' 



IB. 



III. QUANTUM INVERSE SCATTERING FOR SINGLE MODE SPIN-BOSON MODEL 

In this section we apply the general theory presented in the previous section to the spin-boson problem. Starting 
from the i?-matrix (j4]) we consider two types of Lax matrices satisfying ([2|) : the spin degree of freedom is described 
by 

4 S) (A)=f A + ?? t^ " S > I, (17) 



r)Sf X + i] — -qS 



where the operators Sj and Sj 1 are chosen to be irreducible representations of su(2). On way to construct a Lax 
matrix for the bosonic degree of freedom is using a representations of su{2) in terms of Bose operators in (|17p . e.g. the 
Holstein-Primakoff or the Dyson-Maleev representation. Here, however, we choose the bosonic Lax matrix as 



17 



Pi \ 3 \ t , \ i 

If J \ ja-j rja^aj + X + r) 



with [aj,<xt] = Sjk- For the construction of the transfer matrix (|16p we also need the inverse of the bosonic Lax 
operator which reads 



"fdj rjCLjCLj — X + T] 



4 M) (A)=det 9 [4 M) ]-M " ' i \ =det q [£f> 1) }- 1 £f' 2, (-X) . (19) 



Here det g [£^ b,1 ' ) ] = —^(X/rj) is the quantum determinant of C^ ,V> . 
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The i?-matrix {I} intertwines both the spin and bosonic Lax matrices; additionally it is GL(2) invariant, namely it 
commutes with Q®Q for all Q S GL(2) (see ([5])). This GL(2) invariance of the i?-matrix implies that GL(2)-transformed 
Lax matrices G1CG2 are again solutions of the Yang-Baxter relation @. At a formal level, GL(2) transformations 
preserves the structure of the spin Lax matrix (|17p and essentially correspond to a rotation of the local spin. We point 
out that this is not the case for the bosonic Lax matrix (fP8|) ; though general GL(2) transformations lead to canonical 
transformations of the boson degrees of freedom, only diagonal GL(2) transformations preserve the structure of the 
Lax operator. 

In what follows we restrict ourselves to setups with a single spin and a single bosonic mode, corresponding to a 
two-site monodromy matrix. Many relevant features are already contained in this simple case, and an extension to the 
multi-spin and multi-mode situation is technically straight forward. Integrable spin-boson models can be constructed 
e.g. from the transfer matrix for toroidal boundary conditions (periodic boundary with two boundary twists) 

t twis t{\) = tr A [g s c (s \x - z )g b c (b) (\ - (20) 

where Q s ,Qb G GL(2) do not depend on A; the quantities zq,z\, shifting the spectral spectral parameter locally 
are known as 'impurities' in the jargon of QISM. Since [ttwist(X), ttwistifJ-)] — 0, in particular all coefficients of a 
A-expansion of t tW i S t(X) commute with each other. The first order coefficient of this A-expansion of the transfer matrix 
is a linear combination of the spin- and boson operators occurring in the Lax matrix. This operator, where spin- 
and boson degrees of freedom are decoupled, consequently commutes with every Hamiltonian constructable from this 
transfer matrix. This implies that the problem can be effectively reduced to a dressed spin degree of freedom, in which 
the Hamiltonian is block-diagonal [1]. This can be argued from the GL(2) symmetry of the i?-matrix discussed above, 
in that A-independent twists are elements of GL(2) and can hence be absorbed in a redefined spin and boson Lax 
operator. Note that this result is independent of the specific bosonic Lax matrix chosen; either choice of a bosonic 
Lax operator, as given in Eq. (|18j) or obtained from the spin Lax operator (|17|) inserting e.g. the Holstein-Primakoff or 
Dyson-Maleev bosonic representation of su(2), will produce model Hamiltonians with a conserved number operator 
S z + a^a. Consequently all Hamiltonians constructed in this way are of block diagonal form, with more or less 
complicated non linearit y in the bosonic and the interaction part. Examples for such nonlinear generalizations of 



37 



38 



39 



40 



4l| . where both twist matrices were chosen as the same diagonal matrix. 



3- 



H s ^ p h can be found in 

Non-diagonal hermitian twist matrices have been used in 

For obtaining a model including counter rotating terms from the XXX symmetric R matrix, we look at open 
boundary conditions. This amounts to using the transfer matrix (|16p as the generator of integrable Hamiltonians. To 
this end we will consider the following monodromy matrix 

T(A) = K+{\ + r^C^ 1 (-A - z^C^ (-A - Zl )K-(X)C^ (A - zi)£«(A - z Q ) (21) 



where the boundary matrices are taken from equation (j 1 5(1 - 

In the next section we discuss the exact solution of the eigenvalue problem for the transfer matrix t(X) = tr T(A) by 
means of the algebraic Bcthc Ansatz following two approaches. First by 'direct' algebraic Bethe ansatz, i.e. including 
the bosonic Lax operator (|18|) . and then by first solving an auxiliary spin-spin problem and then performing a suitable 
large-spin limit leading to the bosonic degree of freedom. The explicit construction of model Hamiltonians from (j2Tj) 
will be presented in Section IVl 
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IV. ALGEBRAIC BETHE ANSATZ FOR THE SPIN-BOSON TRANSFER MATRIX 



The starting point of the exact diagonalization of the transfer matrix by the algebraic Bethe ansatz is to identify a 
suitable eigenstate of the transfer matrix, the so-called pseudo- vacuum \p) . In the presence of non-diagonal boundary 
matrices, this state can be identified by a method proposed in Refs. 



27 



321 . Sketching their approach, the boundary 



matrices are brought to triangular form by a suitable similarity transformation of the monodromy matrix pip. This 
amounts to finding matrices Mj = 



f (A) = M^T(X)M 

= M^K+iX + ^M^M^C^ 1 (-A - z^M^ 1 ^ 1 (-A - z 1 )M 3 M^ 1 K^(X)M 2 x 
xM 2 " 1 £W(A - zi)Af 1 M 1 - 1 £W(A - z )M 



(22) 
(23) 
(24) 



such that the lower left entry of both Mq 1 K + (X + t])M 5 and M^ 1 K-{\)M 2 vanishes. The remaining parameters in the 
Mj have to be chosen such that the lower left entries of the transformed Lax operators C^ s,h ^ annihilate the spin and 
bosonic vacuum states \lo s ) and \u>b) respectively. As a result we can choose the product state \Q) = \ui s )(gi\u;b) as pseudo 
vacuum for the algebraic Bethe Ansatz. Following Refs. 1271. l32l for the spin Lax matrix we readily obtain x^/yo = xi/yi 



and £4/2/4 = £5/2/5. Note that having M~ 1 £( Q ) Mj-i \u a ) triangular implies the same for Mj-xC^ (A — z )Mj \ui a ), 
where the matrices Mj may differ from Mj in their right column only. Since the transfer matrix is independent of 
this spurious freedom, we are free to choose M3 = M 2 , M4 = Mi, and M5 = Mq without loss of generality. If all Lax 
matrices are spin Lax matrices, all Mj can be chosen to be equal 
spin. 



27 



32j, corresponding to a global rotation of the 



A. Algebraic Bethe ansatz with the bosonic Lax operator 

In this section, we adapt this procedure to the spin-boson model for the two-site monodromy (|21j) and the choice 
£y>) — £( f, > 1 ); the alternative choice — C^ b ^ is discussed below. The approach described next for open boundary 
conditions equally applies to toroidal boundary conditions. 

The lower left entry of the transformed bosonic Lax matrix £<- b \\) = M^ 1 £^{X)M 1 is (up to a rescaling of the 
spectral parameter) 

T t 

- xiy 2 \ + {-x 2 + y 2 a 1 ){r]x 1 a - f3y{) (25) 
r) 

Clearly, this operator can annihilate a vacuum state for all A only if xiy 2 = 0. The bosonic creation operator a) has 
no right eigenstate, therefore we have to choose y 2 = 0. This means that the gauge matrix left to the bosonic Lax 
matrix must already be in upper triangular form. Consequently, the pseudo- vacuum of the bosonic degree of freedom 
is the coherent state {r\x\a — (3yi) \tUb) = with 

£f£ = ££i = £e«\ (26) 



We also comment that for the Eq. (|25|) the structure of two reference states found in 



42| is lost in the spin-boson case. 



Writing C [ > = we have 

AW\u b ) = VL[(\-z 1 )e'>)\u b ) , (27) 

The characteristic factor e p is obtained from the condition that K + = M^ 1 K + M be upper triangular 

e-* = K+ l X (28) 

K+ - #+;22 

where k+ is an eigenvalue of K + . 

Note that the parameters ro, so> ^o> *o as well as r 2 , S2, f2, S2 can be chosen such that both boundary matrices K + and 
K- are diagonal. These parameters enter both the Bethe equations and the expression for the eigenvalue Cambda(X) 
of the transfer matrix, Eq. (|32ll . The transfer matrix and its spectrum, nevertheless, is independent of this choice. 
This freedom provides with a variety of isospectral Hamiltonians and corresponding Bethe equations. The transfer 
matrix, however, differs from that obtained from Eq. (121|) and corresponding diagonalized K + and This is due 
to the fact that for rj, Sj not being all the same, the "gauge" transformations of the Lax operators do not induce 
canonical transformations on the spin- and bosonic degrees of freedom. 

Here, rather than bringing the boundary matrices into diagonal form, we choose gauge transformations to simplify 
the Bethe ansatz equations. To be specific, let detMj = XjSj — yjTj — 1 and Xj = r j = 1 for j = 0, 1 giving 
sj = yj + 1 = c~ p + 1. Furthermore, the condition y 2 = for the existence of the bosonic pseudo vacuum implies that 
the boundary matrix A"_ in (|22p must be upper triangular (when bringing the monodromy matrix in lower triangular 
form instead, the particular role of Xj and yj is taken by rj and Sj respectively). Therefore, we are free to choose 
M 2 = 1 and the action (|27|) of C^(X — z\) on the bosonic coherent state simplifies 

CW|u; 6 ) = 0, 

A^\u b ) = (X- Zl )\co b ) , (29) 
D< 6 > \w b ) = - |w 6 ) . 

Similarly, after the similarity transformation, the action of the spin Lax operator (A — zq) for a 25 + 1 dimensional 
representation of su(2) on \u> s ) becomes 

C^\cu s )=0, 

A^\u s ) = (X-z Q + V (S+l))\co s ) , (30) 
D^\u> s ) = (X-z -t,(S-1))\uj s ) , 
while the boundary matrix K + is transformed into 

-K+, 21 c 2 p + {2K +iU - K+, 22 )cP + 2K+, 12 * 

K+, 21 e 2 P - (K +;11 - 2K+, 22 )eP - 2K+, 12 




± Xy/1 + v + n + * 

£+ T AVI + v+^+ 



(31) 
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It is worth noticing that for a multi-site spin and multi-mode boson generalization of this model the actions of all 
consecutive Lax matrices on the pseudo-vacuum have to be in triangular form for the above procedure to work. 
Furthermore will we briefly discuss the diagonalization procedure, when the bosonic Lax matrix Eq. (|19[) is inserted 
into the transfer matrix, Eq. (fT6|) - and not Eq. fTS]) as above. A straight forward calculation shows that then the 
gauge matrix Mi to the right of the bosonic Lax matrix and all following Lax and boundary matrices must be in 
triangular form. As a consequence, K + instead of K— must be triangular. 

We now proceed to the result of the algebraic Bethe ansatz described above. The eigenvalues A(A;{Aj}) of the 
transfer matrix (|16[) are obtained within the standard formalism [19I ] 



A(A) 



(A + tjS + z )(X + r/{S + 1) - z )] \X - z x 



(A + z Q + nS)(X + z - r)(S + 1)) 
2(A + 77)(A + C-)(A + C + ) 



A + zi 



(2A + V ) 
(X-f 1 S + z Q )(X-f 1 (S- 



J 1 = 1 1 (A-A J + §)(A + A, + f) 



n 



zo) 



(A + z + r)S){\ + z a - r)(S + 1)) 
2A(A-C-+r?)(A-C + +r?) 



A + z\ + 77 



(32) 



(2A + 77) 



X + Z! 

(A-A, + f)(A + A, + f) 



II (A-A, + |) (A + A, + f)' 



The pseudo-vacuum |0) has eigenvalue Ao(A) = A(A; 0). The consequence of choosing C^ b,2 " > as bosonic Lax operator 
consists in an overall minus sign and z\ — > —z\. Eigenstates are characterized by the roots {A^} of the Bethe equations 



(Aj + rjjS - |) + z )(A, + ri(S + go) 

(\i - n(s + i) + z )(x t - rts - 1) - zo) 

(X t - Xj+rfiXi + Xj+r)) 



A, - zi - § 
Xi + zi + f 



(A, + C+ ~ j)(A, +g_ - §) 
(A, -C+ + IXA, -£- + §) 



n 



(33) 



(A; - Xj - 77) (A< + Aj - 77) 



where the boundary matrix i4T + has been rescaled with y/1 + and ^+ := + is the inverse of the 

effective boundary field. Eqs. (|33|) coincide with the Bethe equations for an open chain with diagonal boundary 
matrices after a rotation of the spin and bosonic degrees of freedom. Choosing £ < - b ' 2 ' as the bosonic Lax operator 
induces the change z\ — > — z\. 

We close this paragraph commenting on the general implications of the restriction to triangular K- . Though the 
boundary conditions do break the su(2) symmetry of the system, there is a remaining conserved charge in the 
system in this case. Hence, the spin boson models which can be diagonalized within this approach should not be 
expected to feature counter rotating terms. Nevertheless for non-diagonal triangular the resulting model will 
show some features of counter rotating terms. We emphasize, however, that the triangularity of K- is a constraint 
necessary only for this method for the diagonalization of the transfer matrix which is not related to the integrability 
of the corresponding spin-boson model. 



B. Solution by algebraic contraction of the spin-spin transfer matrix 

It is well known that bosonic degrees of freedom can be realized as spin operators, in a certain limit. Such a limit 
constitutes an example of algebraic contraction [4^]. In this section we demonstrate how this idea can be applied to the 
different stages of QISM method, thus allowing to obtain integrable spin-boson model from certain auxiliary integrable 
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spin model. The approach is pursued here in the case the integrable model is with open boundary conditions, thus 
extending the theory for the toroidal boundaries 22 1 . 



At the level of the Lax matrices, the spin to boson mapping is achieved by the combination of two singular 
transformations, one in the auxiliary and the other in quantum space, to have a finite result in the limit. Explicitly, 
the bosonic Lax matrix is obtained as 22 



C (b) (X) = lim fc(e)£ (J) (A-<5), 



where the matrix fc(e) is a diagonal matrix 



fc(e) = diag 1, ,5 = 7j + V e z /2, J = e^/2 



(34) 



(35) 



and the quantum space transformation is defined by {</~,4yJ + , J 2 } — * ^a* 1 , a, —a* 1 a + 4^|. We have assumed (3 
7 = 77 in the definition of the Bose Lax matrix, Eq. (|18p for convenience. 

We apply the Eq l34l to an auxiliary double-row transfer matrix t a (X) that in the limit of e — > oo is equal to t(X) 



t a (\)=tr\K + (X + v ) [c^(-X-z )\ [£^{-X-6-Z!)\ [k(e)}- 1 K_(X) 

xfc(e)£( J )(A - S - Zl )C (s) (X -z )\. 



(36) 



The above operator t a (X) can be conveniently rotated into i a (X) through a new transformation Uj J7]) acting in 
the quantum spaces, along the line suggested by the property ([6]). This simplifies the the effect of fc(e), acting on the 
auxiliary space. The transformed t a (A) reads 

t,(A) = UUXp- 1 

= tr I k(e)K + (X + r]) (Me))- 1 Lc« (-A - z )] ~* [c^(-A - S - zijl ~* if-(A) 



x C 



(• / )(A-<5-z 1 )£( s )(A-z )} 



(37) 



Therefore, we can obtain the spectral properties of t a (X) by the analysis of t a (X), whose eigenvalues were previously 
obtained by means of algebraic Bethe ansatz In this rephrasing the if-matrix we need to take into account is 



K + (X) = k(e)K + (X)(k(e))-' 
C++ A e 2 Xn+ 



(38) 



Of course, the eigenvalues of the transfer matrix t a (X) and t a (X) are identical and given by the following expression 



A„(A) 



(A + rjS + zq)(X + T) + r/S — zq) 



(r]S + X + zo)(v(S + l)-X-z ) 



(X-5 + rjJ + r]- zi)(A + 5 + rjj + z x ) 



(rjj + X + S + zi){r){J + 1) - A - S - zx) 



2(A + 7 ? )(A + C-)(A + C+) 



(2A + r,) 



n 



(A-A J -f)(A + A,-f) 



i(A-A,+|)(A + A, + |) 



(A - rjS + z )(X + rj - r/S - z ) 
{r)S + X + z )(ri(S + l)-X-z ) 

2X(X-C-+V)^-C++V) 



(2A + 77) 



n 

3=1 



(X-S-rjJ + rj- zi)(A + S - r)J + z x ) 
(r]J + X + S + zi)(rj(J + 1) - A - S - zi) 

j 1 Jil" 1 /s ] 1 "a" 



(39) 



( A -A J + f)(A + A J + f) 
(A-A j + f) (A + A i + 1) 
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while the Bethe ansatz equations are given by 

_ ' (Aj - | + yS + z )(Aj + | + 7?S - z ) 
" (Aj -a-vS + zo){Xj + , j-r 1 S-z a ) 



(Aj + 5 - f + 77J + ^^(Aj- - <J + I + r)J - zx) 



(A, 



6 - § - ryJ + ^(A, - <5 + § - 77J - Zl ) 



(40) 



where C± = 



A,- + C- 



A, + C+ - 



in 



(A 3 - - Aj + t?) (Aj + A, + 77) 
(Xj - Xi - rj) (Aj + A,- - t?) ' 



which do not depend on e. For the algebraic Bethe ansatz, one of the two following 



y/1 + \L±V± 

assumptions on the parameters £±, ji± and v± has been made (see [32[) 



1 + VI + v-\i- 1 + VI + v+v^ 



— V- 



1 + sjl + 



(41) 
(42) 



1 + VI /^-/i- /i+e 
To obtain the eigenvalues and the Bethe ansatz equation of the spin-boson problem, associated to the transfer 



matrix t(X) we take the limit e — > 00 of Eqs. (|4"0"]) . This leads to the same eigenvalue (|32p and Bethe ansatz equations 
(|33f as obtained from the 'direct' ansatz. In the limit of e — > co the constraints (|4T|) . (|42|) become v_ = and = 0, 

respectively. This means that the boundary matrix isT_(A) has to be in a triangular (or even diagonal) form already, 
and £_ = — 

Considering fi± = v± = in the eigenvalue expression and Bethe ansatz equation above we can recover the results 
obtained for the case of diagonal i^±(A) if -matrices without use of the limit process. 

As final remark we observe that a definite limit e — > exists, despite the matrices k{e) are singular in the limit; 
such singularity is compensated by the divergence in the "impurity" 5 (see Eq.(|35jl). 



V. INTEGRABLE HAMILTONIANS 



The transfer matrix t(X) for rational solutions of the YBE is a polynomial expression in the spectral parameter. 
Because the transfer matrix is a commuting family of operators in A, the coefficients of such a polynomial are the 
integrals of the motion of the theory. In the following section we exploit this property of the transfer matrix to extract 
the Hamiltonian. In the section fyBl wc obtain the Hamiltonian following an alternative procedure that is in the spirit 
of the quasi-classical expansion. 



A. Integrable models derived from the transfer matrix 

We analyze the integrals of the motion generated by the transfer matrix (fl"6j) and the bosonic Lax matrix = C^ b ^ 
for spin- 1/2. We therefore expand it in powers of the spectral parameter 

6 

t(A)=:]Tif„A" (43) 

71=0 

where we set rj = 1. 
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There are only two independent constants of the motion in this problem and we find that Hq, H§, and Hq are 
constants and the relations 

2Hi - H 3 = 5n-is + (44) 
H A -H 2 + H 1 = 3/x_i/+ (45) 

With this result for the general case we are now ready for a systematic case study. 

At first, we focus on diagonal boundary matrices K±. In this case we find up to a constant 

Hi = 4/3 7 (S z - n) . (46) 

This constant of the motion readily unveils this Hamiltonian to be of the Jaynes-Cummings type in the RWA because 
the Hubert space splits up in two-dimensional (spin-1/2) subspaces with fixed eigenvalue of Hi. Problems in absence of 
a conserved number operator in general will require more sophisticated techniques as employed for the diagonalization 
of the XYZ chain or the XXZ chain with non-diagonal boundary fields 



34 



For diagonal K + and non-diagonal but triangular _KT_, the expansion (|43|) leads to (up to a constant) 

Hi = Q {S Z - h) + C x ( 1 - + 2S z s ) S + + C 2 a + C 3 a(2S z - h) (47) 



2s 

where s is the spin, Hq = 4/?7, C\ — ^_(2£ + + 2z$ — 1), C 2 = 2^v-{z\ + and C3 = 2jv-. As we will show below, 
the Hilbert space of the non-hermitian Hamiltonian (|47|l is not reduced to two-dimensional subspaces; therefore, it 
describes a resonant spin-boson model (i.e with vanishing de-tuning of the spin-boson frequencies) containing counter- 
rotating terms. It can be treated by the method used above. A careful analysis reveals that the non-hermiticity of the 
Hamiltonian is inherent to this transfer matrix unless the spin degrees of freedom are decoupled from the boson. We 
attribute this feature to the structure of the bosonic Lax matrix (fT5|) used here. 

In what follows we restrict ourselves to spin s = 1/2. A vacuum state for the Bethe ansatz on the Hamiltonian (|47|) 
then is |0, using the notation \n, a) such that S z |n, a) = a |n, a) and h |n, a) = n |n, a). Defining C2 =: C\ + 5, 
C3 = C\ + S + A, the creation operator B(X) has the form 

B(\) = [A 2 \(\ + l)-S 2 + 2SAS z ]{a^(n-^)+h(h-l)} 
-A [25h + 2A(AA(A + 1) - S)a + 2ASah)} S" 

-A(A(A 2 + A + £_)-<5£_)S- (48) 

A right-eigenstate |Ai, . . . ,Xn) '■— -B(Ai) • • •-B(Ajv) |0, 5) involves all states \n,a) with n < N except |JV, —5). It is 
interesting to notice that for positive integer values of the inverse boundary field the Hilbert space accessible from 
the vacuum |0, \ ) is restricted to 

^_eN+ -={\n,a) ; n < ^} (49) 

Without loss of generality we set g\ = 1 and /U+ = v + — 0. The latter corresponds to a diaj 
matrix K+(\) = ^+11 + \a z \ any non-diagonal K+ can be diagonalized by a gauge transformation 27 



onal boundary 
32| and proper 



r 0i so, fo, so- For convenience, we choose 7 = /3 = 1 and rename g 2 =■ g. The most general Hamiltonian generated in 
this way is equivalent to 

H = -^M =: Hs +H b + H sb + H sbb (50) 
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via suitable canonical transformations on the spin and bosonic degrees of freedom. 

Our goal will be to have the spin and boson Hamiltonian together with major part of the spin-boson interaction 
hermitian. To this end, we set v_ — ^—(jj + £+ z i), which makes the occurring terms aS z and a) S z hermitian and the 
choice £_ = 1 — z\ — realizes the same for aS~ and S + . z\ = 1/2 eliminates a term nS + and the pure 

spin Hamiltonian becomes hermitian if 



£+ 2z Q - 6£+ - 1 

The terms a and a' in the pure boson part of the Hamiltonian are obtained hermitian if 



(51) 



Z °=2 ± 



1 Ji-9(g + 2)e + + g 2 (g + i)e + 



(52) 



V9(ge+ - 1) 

and it turns out that both conditions (|51[) and (|52[) are compatible and fixes a one-dimensional manifold in £+-zo 
space (see figure [1]). 



6ft 



- 1 1 1 r- 





Abbildung 1: Left: Solutions of the nonlinear equation in and g emerging from inserting equation (|52|l into (|51[) . Right: the 
corresponding values for zq. 



The resulting Hamiltonian is as follows 

H, = U-^ - 4 1 + 2 ^ 8 e+(2fr-in sz + (€ + -l)((2zo-l) a - 44) 

V 22 + 2£+-l 2z -6e+-iy ' £+(2z -6£ + + l) v 7 

H b = g g£ 2 _ 1 n - -9 g ^ 2 _ 1 (° + a ) + 5 g ^ 2 _ 1 a " ( 54 ) 

ff sb = -8.g£+ nS* + 2^_( 5 C+ + 2) (a + a f )5 z + 2 5 (2z - 2£+ - 1) (aS~ + tfS + ) (55) 
H sbb = gfi- ((2z + 2£+ - l) g ^ 2 a t2 5+ - (2z - 2£+ - 1) a 2 ^ 

+g(2z Q + 2£+ - « 2 5' + + 2 a f n5+) - 4 5 /i_£+ anS* z (56) 



Besides the hermitian parts H s and there are many non-hermitian terms, in particular nonlinear couplings as in 
H a f,b- It is worth focusing on the non-hermitian term in Hf, proportional to an. This term has an immediate physical 
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interpretation: it describes a leakage of bosons/photons from the cavity, where the escape rate is proportional to the 
number of bosons/photons (i.e. the intensity of the cavity boson field). 

B. Integrable models derived from the quasi-classical limit of the transfer matrix 

Another wa y t o obtain an integrable Hamiltonian is to take the so called quasi-classical limit of the transfer 
,l29(. It consists in a series expansion in the 'quantum parameter' r\ of the transfer matrix around r\ = 0: 



matrix 



20 



21 



t(A) = f(°) +t]t( 1 ' (\) + r] 2 (A) + . . . with the aim of creating a commuting family of quasi-classical transfer matrices 
7"' fe )(A). This procedure may be particularly useful for extracting 'simple' though non-locally interacting Hamiltonians 
out of the transfer matrix. Examples are the Gaudin magnets and corresponding BCS-like models. 

There is a wide freedom of introducing an 77-dependence to the boundary matrix parameters in the theory, without 
destroying its integrability. For the case of integrable spin chains with diagonal boundaries [2^] this procedure was a 
valid tool to define integrable one-parameter extension of Gaudin models in non uniform magnetic fields 25]; such 
models describe metallic grains with pairing and magnetic interaction 24j . 

In the present case the transfer matrix is the finite sum 

t(X) = f 1 - k f^-^ + ■■■ + f ( °) + 77 f « (A) + t] 2 f (2 > (A) + • • • + r, m T (m) (A) (57) 
where k and m are integers. Expanding the commutator relation [t(X),t(X')] = in 77, we obtain 

2m 

[*(A),t(A')]= £ r 1 l C l (X,X') = 0, 

l = ~2k 

which implies C;(A, A') = for all I, The first relevant terms are 

C_ 2fc (A,A')=[f(- fc )(A),f(- fc )(A')], 
C_ 2fc+1 (A,A') =[f(- fc )(A),f(- fc+1 )(A')] + [f(- fe+1 )(A),f(- fe )(A')] , 

C- 2fe+2 (A,A') =[f(- fc )(A),f(- fc+2 )(A')] + [f(- fe+2 >(A),f(- fc )(A')] + [f(- fc+1 )(A),f(- fc+1 )(A')] , 

(58) 

^From the expressions above, one finds that the first f^ n \X) which is not a C- number (times the identity) gives 
rise to a family of commuting operators. Generically, the lowest order f (~ fc ) is a C-number. Therefore the first class 
of integrable models were generated by [f( -fc+1 ) (A), f(~ fc+1 ) (A')] = 0. In the presence of boundary matrices, these 
are typically non-trivial operators but representing non-interacting Hamiltonians. The task is then to tune the free 
parameters such that f( _ ' t ' +1 '(A) is also a C-number, and that the lowest non-trivial order in 77, e.g. f(~ fe+2 )(A), be an 
interesting Hamiltonian. 

Here we allow £±, p,±, v± to be generic functions of rj whose form can be fixed to modify to our convenience the 
^-expansion of the transfer matrix, and ultimately the Hamiltonian. Specifically, we write the boundary parameters 
as 

x± = a^V 1 + ^± 0) + x± )r n (59) 

with x = {£,/!, k}, and set ^ = 2/£ 1 ^ = £.1 ^ = 0. Then the lowest non-vanishing term in the expansion, 
proportional to r/ -2 , already contains non trivial operators; they can be made vanishing by the choice fffl = v 1 ^ = 0. 
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Analogously, the next term in the expansion can be made proportional to the identity by = $P = and z\ = 0. 
The next non-trivial term in the 77-expansion (A) can be used as Hamiltonian 

H = t^t (0) (A) = 4Q <x f a + 2a (a + a*) + A Z S Z + 4jz aS x 

+GS 2 a t a - 4720(05"- + a^S+) + 7 2 + S~ S + - 2{S z f] + Ct (60) 

where il = -(\ 2 -z%), a = {3v ( ^ - 7 i/| 1} and A z = -4A 2 - 2 7 2 (1 - 2^ f 1) ), G = 4 7 2 (1 - 7) and C = 
- 1Zl {n [2(2 + C- ~ «M)A 2 + 2£° A+] - 2(A - z )( 7 - A^ )}- 
For S = 1/2 the Hamiltonian simplifies to 

fli/a - ^Wa + A,S* + GS z a)a + GaS z (a + a f ) - .g (5+^ + S~a) (61) 

where the bosonic operators have been displaced a —>■ a + a/2, A z = A z + a 2 G and a term proportional to the identity 
was omitted. The Hamiltonian (|61|) display the rotating and counter-rotating simultaneously. This can be evidenced 
by a rotation R of the spin 



RH lj2 R- x = Q a^a + A z [cos(26)S z - sm(29)S x ] + G [cos(26)S z - sm(29)S x ] a) a 

cos(26>) + 1 . , + n _ s cos(26») - 1 ^ n 

-5 ~y {S + a'+S a) - g (S + a + S a 1 ) . 

where R = exp[6(S + — S~)] and tan^ = Ga/g has been chosen to eliminate the S z (a + <r) term. 



VI. CONCLUSIONS 



By the Quantum Inverse scattering method we have constructed integrable models with twisted and open boundary 
conditions which contain rotating as well as counter rotating terms. 

Twisted boundary conditions lead to Hamiltonians where solely rotating (or counter-rotating) terms appear. 

In case of open boundaries, we diagonalize the models by algebraic Bethe ansatz when at least one of the boundary 
matrices has triangular form. To this end, we extended the procedure for spin chains to models which are derived from 
a particular bosonic Lax matrix besides the standard spin Lax matrix. We further analyzed generic boundary matrices 
in the sense that they cannot be brought both in triangular form at the same time. For this choice of boundary fields, 
the symmetry leading to conservation of the 'number operator' S z + ^ . aja^ is broken. 

Taking the full transfer matrix the integrable Hamiltonians are non-hermitian. However, a parameter choice is 
found that shifts the non-hermiticity to non-linear terms in the spin-boson interaction and the bosonic part of the 
Hamiltonian. In the latter, the non-hermiticity can be interpreted as a density dependent leakage of photons out of 
the cavity. Different parameter choices could be analyzed in order that non-hermitian terms will appear elsewhere, 
e.g. in the spin part of the Hamiltonian. It is worth studying what type of non-hermiticity are most realistic in order 
to tailor the parameters closest to this situation. The diagonalization of the resulting Hamiltonians together with a 
discussion of the effects of the undesired non-hermitian parts is left to forthcoming work. 

Interestingly enough, by a procedure that is close in spirit to the quasi-classical limit, also hermitian Hamiltonians 
can be extracted out of the transfer matrix. We have given an explicit example, which contains both rotating and 
counter-rotating terms. The evidence we have for this is that any eigenstate (i.e. coherent states) of the bosonic part 
is spread out all over the bosonic Hilbert space by means of the interaction term. 



16 



Though integrable counter-rotating spin-boson models have been formulated in this work, it must be emphasized 
that their diagonalization with e.g. the algebraic Bethe ansatz is still an open problem. It could be interesting to 
study general open boundary conditions still for the XXX model but based on the Holstein-Primakoff representation 
of su(2). The advantage of such an approach is that the Hamiltonian is constructed in a straight forward manner from 
the existing solution of the open XXX spin chain Hamiltonian and hence also the resulting spin-boson Hamiltonian is 
hermitian by construction. The challenge again lies in the diagonalization of the transfer matrix. Finally, the analysis 
of spin-boson models fanning out from a bulk with XYZ symmetry (parametrized by elliptic functions) constitute 
an interesting challenge for future investigation. 
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